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Abstract 

A fully canonical quantization of electromagnetic field is intro- 
duced in the presence of an anisotropic polarizable and magnetizable 
medium . Two tensor fields which couple the electromagnetic field 
with the medium and have an important role in this quantization 
method are introduced. The electric and magnetic polarization fields 
of the medium naturally are concluded in terms of the coupling tensors 
and the dynamical variables modeling the magnetodielectric medium. 
In Heisenberg picture, the constitutive equations of the medium to- 
gether with the Maxwell laws are obtained as the equations of motion 
of the total system and the susceptibility tensors of the medium are 
calculated in terms of the coupling tensors. Following a perturbation 
method the Green function related to the total system is found and 
the time dependence of electromagnetic field operators is derived. 
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1 introduction 



One of the most important quantum dissipative systems is the quantized 
electromagnetic field in the presence of an absorbing polarizable medium. In 
this case there are mainly two quantization approaches, the phenomenological 
method [l]-[7] and the damped polarization modeljElE]. The phenomenolog- 
ical scheme has been formulated on the basis of the fluctuation- dissipation 
theorem [10]. In this method by adding a fluctuating noise term, that is the 
noise polarization field, to the classical constitutive equation of the medium, 
this equation is taken as the definition of the electric polarization operator. 
Combination of the Maxwell equations and the constitutive equation in the 
frequency domain, gives the electromagnetic field operators in terms of the 
noise polarization field and the classical Green tensor. A set of bosonic op- 
erators is associated with the noise polarization which their commutation 
relations are given in agreement with the fluctuation- dissipation theorem. 
This quantization scheme has been quite successful in describing some elec- 
tromagnetic phenomena in the presence of a lossy dielectric medium [TTj-|15]. 
The phenomenological approach has been extended to a lossy magnetic or 
anisotropic medium [16] -[18]. This formalism has also been generalized to an 
arbitrary linearly responding medium based on a spatially nonlocal conduc- 
tivity tensor [T9] . 

The damped polarization model to quantize electromagnetic field in a disper- 
sive dielectric medium [8j [9] is a canonical quantization in which the electric 
polarization field of the medium applies in the Lagrangian of the total system 
as a part of the degrees of freedom of the medium. The other parts of the 
degrees of freedom of the absorbing medium are related to the dynamical 
variables of a heat bath describing the absorptivity feature of the medium. 
In this method the dielectric function of the medium is found in terms of 
the coupling function of the heat bath and the polarization field, so that it 
satisfies the Kramers-Kronig relations [20]. This quantization method has 
been generalized to an inhomogeneous medium [2~Tj. 

In the present work we generalize our previous mo del [2"2"l |2"5] to an anisotropic 
dispersive magnetodielectric medium using a canonical approach. In this 
formalism the medium is modeled with two independent collections of vec- 
tor fields. These collections solely constitute the degrees of freedom of the 
medium and it is not needed the electric and magnetic polarization fields to 
be included in the Lagrangian of the total system as a part of the degrees 
of freedom of the medium as in the Huttner-Barnett model [HJ. In fact the 
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dynamical fields modeling the dispersive medium are able to describe both 
polarizability and absorptivity features of the medium. 
This paper is organized as follows. In Sec. 2, a Lagrangian for the to- 
tal system is proposed and classical electrodynamics in the presence of an 
anisotropic polarizable and magnetizable medium with spatial-temporal dis- 
persion briefly is discussed . In Sec. 3, applying the Lagrangian introduced 
in Sec. 2 a fully canonical quantization of both electromagnetic field and 
the dynamical variable modeling the responding medium is demonstrated. 
Then in Sec. 4, the constitutive equations of the medium are obtained as the 
consequences of the Heisenberg equations of the total system and the electric 
and magnetic susceptibility tensors of the medium are calculated in terms of 
the parameters applied in the theory . In Sec. 5, it is shown that the Green 
function of the total system in reciprocal space satisfies an algebraic equa- 
tion and a perturbation method to obtain the Green function is introduced. 
Finally in section 6 the model is modified for media which the distinction 
between polarization and magnetization is not possible. This paper is closed 
with a summary and some concluding remarks in Sec. 7 . 



2 Classical electrodynamics in an anisotropic 
magnetodielectric medium 

In order to present a fully canonical quantization of electromagnetic field 
in the presence of an anisotropic polarizable and magnetizable medium, we 
model the medium by two independent reservoirs. Each reservoir contains a 
continium of three dimensional harmonic oscillators labeled with a continu- 
ous parameter u. We call these two continuous sets of oscillators " E field " 
and " M field" . The E field and M field describe polarizability and magne- 
tizability of the medium, respectively. This means that, in this approach it 
is not needed the electric and magnetic polarization fields of the medium to 
be appeared explicitly in the Lagrangian of the total system as a part of the 
degrees of freedom of the medium , but the contribution of the medium in the 
Lagrangian of the total system is related only to the Lagrangian of the E and 
M fields and these fields completely describe the degrees of freedom of the 
medium. The presence of the "E field" in the total Lagrangian is sufficient 
for a complete description of both polarizability and the absorption of the 
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medium due to its electrically dispersive property. Also the " M field" solely 
is sufficient in order to description of both magnetizability and the absorp- 
tion of the medium due to its magnetically dispersive property. Therefore, 
in order to have a classical treatment of electrodynamics in a magnetodielec- 
tric medium, we start with a Lagrangian for the total system (medium + 
electromagnetic field ) which is the sum of three parts 



L(t) — L res + L em + L 



int 



(1) 



where L res is the part related to the degrees of freedom of the medium and 
is the sum of the Lagrangians of the E and M fields 



where 



and 



du I d 3 r 



L e + L r 



-A^ • X w 



-oj 2 X u ■ X u 



du / d 3 r 



_ - 1 fil J in 



u 2 Y„, ■ Y„ 



(2) 
(3) 
(4) 



Here the fields X u and Y^ are the dynamical variables of the E and M fields, 
respectively . 

In (CQ) L em is the contribution of the electromagnetic field in the La- 
grangian of the total system 



err 



-e E 2 

2 2/i 



(5) 



and L int is the part describing the interaction of the electromagnetic field 
with the medium 



J int 



du / d 6 r / d 6 r'f ij (u,r,r')E\r,t)X 3 ul (r',t) + 



du I d A r / rf 3 r / £ / ij (o;,f,r / ) J B i (f,t)yj(r , ,t) 



(6) 



The contributions L e and L m in the Lagrangian L res are equivalent to the 
consequences of diagonalization processes of the matter fields in the Huttner- 
Barnet model [8]. That is, L e ( L m ) is equivalent to the diagonalization 
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of the contributions of related to three parts in the Huttner-Barnet model: 
the dynamical variable describing the electric polarization (magnetic polar- 
ization)of the medium , a heat-bath B ( B') interacting with the electric 
polarization (magnetic polarization ) and the interaction term between the 
heat-bath B ( B') and the electric polarization (magnetic polarization). In 
the present approach modeling the medium , in a phenomenological way, 
with two independent set of oscillators the lengthy diagonalization processes 
have been eliminated in the start of this quantization scheme. Particularly 
the diagonalization processes may be more tremendous for an anisotropic 
medium. 

In equations flS]) and flS]) E = — ^ — V<£> and B = V x A are electric 
and magnetic fields respectively, where A and ip are the vector and the scalar 
potentials. The tensors / and g in (jSJ) are called the coupling tensors of the 
medium with electromagnetic field and for an inhomogeneous medium are 
dependent on the both position vectors r and r'. The coupling tensors are 
the key parameters of this theory. As was mentioned above, it is not needed 
the electric and magnetic polarization fields of the medium explicitly to be 
appeared in the total Lagrangian ([I])-© as a part of degrees of freedom of 
the medium . As we will see, the electric polarization( magnetic polarization 
) of the medium is obtained in terms of the coupling tensor / ( g ) and the 
dynamical variables X w ( Y u ). Also the electric susceptibility tensor ( mag- 
netic susceptibility tensor) of the medium naturally will expressed in terms 
of the coupling tensor / ( g ). The coupling tensors / and g are appeared as 
common factors in both the noise polarization fields and the susceptibility 
tensors of the medium, so that for the free space the susceptibility tensors 
together with the noise polarizations become identically zero and this quanti- 
zation scheme is reduced to the usual quantization of electromagnetic field in 
free space. Furthermore when the medium tends to a non-absorbing one, the 
coupling tensors and the noise polarizations tend also to zero and this quan- 
tization method is reduced to the quantization in a non-absorbing medium 
[221121]. 

In order to prevent some difficulties with a non-local Lagrangian such 
as in it is the easiest way to work in the reciprocal space and write all 
the fields and the coupling tensors / , g in terms of their spatial Fourier 
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transforms. For example the dynamical variable X u can be written as 

1 



X u {r,t) 



d 6 kX^(k } t) e lk ' r 



(7) 



Since we are concerned with real valued fields in the total Lagrangian ([I])-©, 
we have X L Jk, t) = X^(—k, t) for the field X u (f, t) and the other dynamical 
fields in this Lagrangian. Similarly the real valued coupling tensors / and g 
can be expressed in reciprocal space as 

1 



fij(u),r,r') 



(2tt) 3 
1 



d 3 k I d 3 k'l i .(u J k,k')e lkr '- lk ' r ' 



d 6 k / d 3 k'g..(u,k,k')e 



*/\ ik-f—ik'-r' 



(2tt) 3 

which obey the following conditions 

g (u,k,k f ) = g* (u,-k,-W) 



(8) 



(9) 



■-ij —tj 

The number of independent variables can be recovered by restricting the 
integrations to the half space k z > 0. The total Lagrangian (CO)-® is then 
obtained as 

W) = L res (t) + L em (t) + L mt (t) (10) 



du / d^hX^l 2 -u 2 \X^\ 2 )+ du / d 3 k [\Y u \ 2 -u; 2 \Y 



oo /■/ 



Len.it) = / d 3 k [ E \A[ 2 + S \k if\ 2 - 



\k x A\' 



(11) 

+e I d 3 k ( -ik ■ A_L(f + h.c 

(12) 



L int (t) 



du J d 3 q I d 3 p [(A{q, t) + iq <p{q, t)j ■ f(u, -q, p) ■ X^{p, t) + h.c 
du [ d 3 q [ d 3 p\(jt(q,t)-tq^(q,t)) ■ f(u,q,p)-X^{p,t) + h.c 



+ / du d 3 q d 3 p 



iq x A(q,t)) -g(u,-q,p) -Y^(p,t) + h.c 



du / d 3 q / d 3 p l-iqx A*(q,t)) ■ g(u,q,p) -Y^(p,t) + h.c (13) 



where J d 3 k implies the integration over the half space k z > ( hereafter , 
we apply the symbol / d 3 k for the integration on the half space k z > and 



d 3 k for the integration on the total reciprocal space). In the reciprocal 

space the total Lagrangian ffTO]) - ffT3|) do not involve the space derivatives of 
the dynamical variables of the system and the classical equations of the mo- 
tion of the system can be obtained using the principle of the Hamilton's least 

action, 8 I dt L(t) = 0. These equations are the Euler-Lagrange equations. 



For the vector potential A(k,t) and the scalar potential tp(k,t) we find 

5L \ 5L 

i = 1,2,3 



5{£(k,t))J 8(A*(k,t)] 
fjL £oA(k, t) + /i £o ik (p(k, t) — k X Ik x A(k, t) 
Ho P(k,t) + ifi k x M_(k,t) (14) 



8L 8L 
5(g(k,t)) ) S{ip*(k,t)^j 

=^ -e ik ■ A(k,t) + £ \k\ 2 (p(k,t) = -ik ■ P{k,t) (15) 
for any wave vector k in the half space k z > where 

du J d 3 pl(u,k,p)- X^{p,t) (16) 

M(k,t) = du J d 3 p g(u,k,p) ■ Y^(p,t) (17) 

are respectively the spatial Fourier transforms of the electric and magnetic 
polarization densities of the medium and it has been used from the relations 
(Q. Therefore in this method the polarization fields of the medium are 
naturally concluded in terms of the coupling tensors /, g and the dynamical 
variables of the E and M fields. Similarly the Euler-Lagrange equations for 
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the fields and for any vector k in the half space k z > are easily 
obtained as 

d ( 5L \ SL 

— I —. = v \ = ^ = t=l,2,3 

2L(k, t) + uj 2 X^(k, t) = - J d 3 q f(u, q, k) ■ [X{q, t) + iq^q, t)) 

(18) 



— I 5 ~ I S ~ - i - 1 2 3 

=>- £«(M) £«(M) = J d 3 q g j (co,q,k) ■ (i qx A(q,tj^ 

(19) 

where f^andg^ are the hermitian conjugate of the tensors / and g, respec- 
tively. 



3 Canonical quantization 

Following the standard approach, we choose the Coulomb gauge k-A(k, t) = 
to quantize electromagnetic field. In this gauge the vector potential A is a 
purely transverse field and can be decomposed along the unit polarization 
vectors e A £ A = 1 , 2 which are orthogonal to each other and to the wave 

— * 

vector k. 

2 

A(k,t) = Yl Mk,t)e xU (20) 

A=l 

Although the vector potential is purely transverse, but the dynamical fields 

— * — * 

X w and may have both transverse and longitudinal parts and can be 
expanded along the three mutually orthogonal unit vectors e A £ A = 1,2 
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and e 3 £ = k — t|t as 



A=l 
3 

to(k,t) = Y^^x(k,t)e^ (21) 



A=l 



Furthermore in Coulomb gauge the Euler- Lagrange equation (j!5p can be 
used to eliminate the extra degree of freedom ip(k,t) from the Larangian of 
the system 

e \k\ 2 

The total Lagrangian (II 01) - (II 31) can now be rewritten in terms of the indepen- 
dent dynamical variables A x A = 1,2 and J lilA , Y^ A A = 1,2,3 which 
constitute completely the coordinates of the total system 



d 6 k 



\ik-P\ 2 



+ / d 6 k 



\k\ 2 

EMfe ) • + x E^e A H • K + h.c 



A=l / \ A=l 



(23) 



where the polarizations P_ and M have been defined previously in terms of the 
dynamical variables of the E and M fields in equations (Tl6|) and (fl7|) . The 
Lagrangian ( 12 3 p can now be used to define the canonical conjugate momenta 
of the system. For any wave vector k in half space k z > this momenta are 
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defined as 



~D x (k,t) 

ILa(M) 



SL 



S[A x (k,t) 
SL 



5(2La(M) 
SL 



e A x {k,t) -e Ajr P(k,t) 
= 2Lx(k,t) 



(24) 



The total system are quantized canonically in a standard method by imposing 
equal-time commutation relations between the coordinates of the system and 
their conjugates variables as follows 



Al(k,t) , -D x ,{k>,t) 



ihS xx >S(k — k') 

lhS XX 'S(u) — L)')6(k — k') 

ihS xx >S(cu — co')S(k — k') 



(25) 



Using the Lagrangian (|23|) and the conjugates momenta introduced in (j2 
the Hamiltonian of the total system can be written in the form 



H(t) 



d A k- =- 



e Q \k\ 



tkxY^AxHk) -K+h.c 

\ A=l / 

du j d^{\Q^ x \ 2 + ^ 2 \^x\ 2 + \TLx\ 2 + ^\^ 

A— 1 



(26) 



10 



3.1 Maxwell equations 

Using Heisenberg equations for the operators D_ x and A x , that is 

R x (k,t)-e xl : ■ P(k,t) 



A x (k,t) = -\H ,Mk,t) 



CO 



R x {k,t) 



h 



H ,D x {k,t) 



\k\ 



-A x (k,t)-e 



ik x M(k,t) 



(27) 



Maxwell equations can be obtained as is expected. Multiplying both sides of 
Eqs fl27|) in the polarization unit vector e A £ and then summation on polar- 
ization index A = 1, 2 and using fl22|) clearly give us 



D(k, t) = e E(k, t) + P(k, t) (28) 

3(k,t) = ikxH (29) 

2 

where D = D->fixk pl a Y s the role of the displacement field, E_ = —A — 

A=l 

k{k-P) . 



zk(p 



-A 



so 



is the total electric field and H_(k, t) = — — M is 



the magnetic induction field. 



/'o 



4 The constitutive equations of the medium 

It is well known that the constitutive equations of a responding medium are 
the consequences of the interaction of the medium with electromagnetic field. 
Any quantization method such as this theory, in which the medium enter di- 
rectly in the process of quantization and the interaction of the medium with 
the vacuum field is explicitly given, must be able to give the constitutive 
equations of the medium using the Heisenberg equations of the total system. 
Also the susceptibility tensors of the medium which are a measure for the 
polarizability of the medium should be specified in terms of the parameters 
describing the interaction of the medium with electromagnetic field. In this 
section, we find the constitutive equations of the medium using the Heisen- 
berg equations of the dynamical fields modeling the medium. As we do this, 
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the electric and magnetic susceptibility tensors of the medium are naturally 
found in terms of the coupling tensors / and g. In the Heisenberg picture, 
using the commutation relations fl25|) and the total Hamiltonian fl26|) . the 
equations of motion for the canonical variables X^x , Q , A = 1,2,3 for 
each wave vector in the half space k z > follow as 

X,a(M) = ±[H,X^$,t)) = QJ$,t) 

Q^ x (k,t) = [H,Q^ jX (k,t)] = -u J 2 X^ x (k,t)+ J d 3 qf..(u,q,k)E i (q,t)e j ^ 

(30) 

u qU-l(q,t)) 

where E_(q, t) = —A(q, t) is the Fourier transform of the 

total electric field. Multiplying both sides of Eqs. (130]) in the mutually 
orthogonal unit vectors e A £ and then summation on A = 1, 2, 3 and using the 

3 

completeness relations e Afc e Afc = ^ we ^ nc ^ 

A=l 

X(k, t) + u 2 X(k, t) = J d 3 q P(u,q,k) ■ E{q,t) (31) 
This equation can be integrated formally as 

X^(k,t) = Q (k,0)^^- + X^(k,0)cosut 

+ f dt ^^-V f d s q f 1 {Uj % k) ■ E(q, (32) 
Jo u J 

Now, by substituting X L J k, t) from (1321 in the definition of electric polariza- 
tion density given by ( TT6|) . we find the constitutive equation of the medium 
in reciprocal space relating the electric polarization to the electric field 

P(k, t) = P N (k, t) + e J ' dt' j d 3 q x e (k, q, \t\ - f) ■ E(q, ±t') (33) 

where the upper (lower) sign corresponds to t > (t < 0). Here the \ e ls 
the spatial Fourier transform of the electric susceptibility tensor 

±J™chj^fd s pf(uXp)-f(oJ,q,p) t>0 
X e (k,q,t) = { (34) 

t < 



12 



and is guaranteed to possess solutions for the coupling tensor / in terms of 
the temporal Fourier transform of x e using a type of eigenvalue problem [19j. 
It must be pointed out that for a given x e the coupling tensor / satisfying the 
relation (jMj) is not unique. In fact if / satisfy fl34l) . for a given susceptibility 
tensor ,then 

f'(u,k,q) = I d 3 p' f(oo,k,p') ■ AMJ) (35) 



also satisfy the ( ]3"%|) where A is a tensor with orthogonality condition 

J d 3 pA{u,p,p') ■ A}(u,p,p") = I5 3 (p' -p 1 ') (36) 

Although for a given x e various choices of the tensor / satisfying ( [33]) af- 
fect the space-time dependence of electromagnetic field operators, but all of 
these choices are equivalent and the commutation relations between electro- 
magnetic field operators remain unchanged [23]. That is, the commutation 
relations between electromagnetic operators and the physical observables fi- 
nally are dependent only on the given susceptibility tensor and not on the 
coupling functions f,g. 

In (13"3"j) . P_ N is the spatial Fourier transform of the noise electric polarization 
density which is expressed in terms of the the dynamical variables of the " E 
field" at t = 

P N (k, t) = jf du J d 3 p l(co, k,p) ■ (q^ip, 0)^-^ + XM 0) cosu?j 

(37) 

In a similar fashion the constitutive equation relating the magnetic polariza- 
tion density of the medium to the magnetic field B_(q, t) = iqx A(q, t) can be 
obtained straightforwardly using the Heisenberg equations for the conjugate 
dynamical variables Y^ x and 11^ as the form 

M(k, t) = M N (k, t) + — /' ' dtf [ d 3 q X m (k, q, \t\ - t') ■ B(q, ±t r ) (38) 

where the x m i s the spatial Fourier transform of the magnetic susceptibility 
tensor which is written in terms the coupling tensor g as 

X m (k,q,t)={ ' (39) 

t < 
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and 



M N (k,t) = jf duj j d 3 pg(u,k,p)- (^(p, 0)^p + Y^(p, 0) coswij 

(40) 

is the noise magnetic polarization density. Therefore modeling an anisotropic 
magnetodielectric medium with two independent collections of harmonic os- 
cillators, that is, the E and M field, we successfully have constructed a 
fully canonical quantization of electromagnetic field in the presence of such a 
medium, so that the Maxwell equations together with the constitutive equa- 
tions of the medium have been obtained from the Heisenberg equations of 
the total system. In this method the susceptibility tenors of the medium are 
concluded in terms of the coupling tensors which are the key parameters of 
this theory and are describing the coupling of the electromagnetic field with 
the medium. 



The Equations ([54"]) and (159]) relate the parameters of this theory, that is 
the coupling tensors / and g, to the physical quantities \ e an d X m ■ If the 
coupling tensors / and g are specified, so that the right hands of (154"]) and 
(159]) become identical to the electric and magnetic susceptibility tensors of 
the medium, then according to the constitutive equations (1551) and (1551) the 
operators P_ and M_ defined by ( fl6l) and ( fTTl) are respectively the polarization 
and magnetization of the medium. 

There are media for which the polarization ( the magnetization) is depen- 
dent on both the electric and magnetic field. This quantization scheme can 
be generalized for such media if in the interaction Lagrangian ([6]) each of 
the dynamical variables X u and Y u is interacted with both the electric and 
magnetic field. 



5 Space -time dependence of the electromag- 
netic field operators 

Our goal for this section is to determine the complete expressions of the 
electromagnetic field operators for both negative and positive times. The 
Heisenberg equations of the total system, that is the Maxwell equation ( ]29j) 
together with the constitutive equations ( 128]) . (155j) and (158j) . constitute a 
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set of coupled linear equations which can be solved in terms of the initial 
conditions by the temporal Laplace transformation. This technique has been 
used previously in the damped polarization model[21J. For an arbitrary time- 
dependent operator T(t) the backward and forward Laplace transformations, 
denoted respectively, by T b (s) and T^(p) are defined as 



rf( s ) 



dt e~ st T(-t) 



dt e~ st T(t) 



(41) 



Carrying out the backward and forward Laplace transformations of Eqs. (l28l 

. (j251). (133)1 . (1551) and ik x E(k,t) = -R(k,t) and then their combination 
for the electric field operator we find 



d q k x p,(k, q,s)[qxE_(q,s 



-b 



d q e(k, q, s)E (q, s) = J_ N (k, s 



d 3 q k x fi(k, q,s) (qx E (q, s)) - — / d 3 q e(k, q, s)E (q, s) = J N {k, s) 

(42) 



where 



i{k,q,s) = S(k-q)I + x e (k,q,s) 
fi(k,q,s) = 5(k-q)I-x m (k,q,s) 



(43) 



are respectively the permitivity and permeability tensors of the medium in 
Laplace language. The source terms in the right hand of the inhomogeneous 
Eqs.( l42p are the backward and forward Laplace transformations of the noise 
current density and are expressed in terms of initial conditions of the dynam- 
ical variables of the total system at t = as 



J_^(k, s) 



J_jy(k, s ) 



Hos 2 E N (k, s) - fx s ik x M b N (k, s) + 
ikx d 3 q jl(k, q, s) B(q, 0) — ^sD_{k, 0) 



~~*f — * > t 

I^qs 2 E-N^k, s) + fj,QS ik x M N (k, s) 

-ik x / d 3 q jl(k, q, s) B(q, 0) - fi sD(k, 0) (44) 
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(45) 



with 



A {j (k, q, s) = eis^SafSj^q 13 ^ 01 ^, q, s) 



c 



^ij{k,q,s) 



(46) 



From (l4"5l) we see that the backward and forward Laplace transformations of 
the electric field operator satisfy an inhomogeneous equation with a source 
term. To solve such equations it is common to use the Green function 
method. The suitable Green function associated to the Eq. ( 1451) is defined 
as the solution of an algebraic equation as follows 



Accordingly the solution of Eqs. ( 1431) in terms of the Green function defined 
in (14T|) can be written as 



In the special case of a homogeneous medium for which the electric and 
magnetic susceptibility tensors are translationally invariant, that is when the 
tensors x e {^i r 'it) an d X m (^*! r '?^) are a function of the difference r — r' and 
so is then, the permitivity and permeability tensors e(r,r',t) and fi(r,r',t), 
we deduce 




(47) 




(48) 




(49) 



and therefore from the Eqs. (I46p and fl4T|) we find 
G{k,p,s) =T-\k,s)5(k-p) 




(50) 



16 



where it should be summed on the indices a, f3, 7, 5 in the first term in the 
bracket. For a general medium it may be impossible to solve Eq. (j47p exactly 
to obtain the Green function G. However one can use a suitable iteration 
method and find the Green function up to an arbitrary accuracy which may 
be useful for media with considerably weak polarizability. Let us write the 
tensor A(k, q, s) as the sum of two parts 

A(k, q, s) = A<°>(fc, q, s) + A«(£, q, s) (51) 

where 



A_g } (k, q, s) = ( hkj - 5 l3 k 2 - S -5 t] ) S(k -q) = L(k, s)S(k - q) (52) 



9 2 

6 

is the tensor A in free space, that is when there is no medium and 



s 2 

Aij (k, q, s) = -e iSl e al3j ksqpX™(k, q, s) - — x* (k, q, s) (53) 

is the part due to the presence of the medium, that is the effect of the 
medium in the total tensor A. Then, following the Born approximation 
method previously used in the scattering theory [25], the Green function G 
can be expressed as a series as follows 

G(k,p, s) = Gj°\k,p, s) + G^(k ,p,s) + G ( - 2 \k,p,s)+ ■■■ (54) 

where 

G}°\k,p,s) = L-\k,s) S(k-p) (55) 

is the Green function for the free space and G {n) (k,p,s) n = 1, 2, 3 ■ ■ • are 
found from the following recurrence relation 

d 3 q A<®(k,q,s) ■ Q} n \q,p,s) = - J Sq A^\k,q,s) ■ G} n ~ l) (q, p, s) 

G (n) (fc, p, s) = -L^ (jfe, s) [ d 3 q A (1) (jfe, q, s) ■ Q} n ~ l) (q, p, s) (56) 



where we have used Eq. (l52l . Using the series ( I54jl and the recurrence relation 
( 156]) one can obtain the Green function G up to an arbitrary accuracy. In 
n'th order approximation it can be neglected from the terms j > n + 1 

in the series (|54|) which may be useful for media with considerably weak 
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polarizability. 

Now, having the Green function G, the time dependence of electric field 
operator for negative and positive times can be obtained by inverse Laplace 
transformations of E_ and i£ , respectively. Carrying out the inverse Laplace 



transformation of Eq. (jlSj) for t < we deduce 

1 r+°° r t 

E(k,t) = — du e~ wt \ d 3 p G(k, p, iu + 0+) • J N (p, iu + 0+) (57) 
2vr > /_ 00 J 

where + is an arbitrarily small positive number. For t > the inverse 
Laplace transformation of E_ is 

- - 1 f + °° f -4 

E(k,t) = — / due~ luJt / d 3 pG(k, p, -ioj + + )- J N (p, -iu + + ) (58) 

Because the Laplace transformations of the susceptibility tensors are analytic 
functions respect to the variable s in any point of the half-plane Re[s] > 
then, clearly the cu-dependence integrand in Eq. ( IBTj) also is analytic in the 
half-plane Im[u] < 0. Accordingly, the integral over uj in (157j) vanishes for 
the positive times . Likewise the integral in f l58|) is zero for t < 0. Therefore 
we can combine the two expressions fl57|) and fl58|) into a single one and write 
the full time-dependence of electric field for all times as 



1 



— oo 



E(k,t) = — I du e~ wi I d 3 p G(k, p, %u + + ) • J N (p, iu + 0+) 



-6 



+ G(k, p, -iu + + ) ■ f N (p, -%u + + )j (59) 

It is remarkable that using the asymptotic behavior of the Green function G 
and the susceptibility tensors for large \u\it can be verified straightforwardly, 
the integral in (159]) reduces to E_(k, 0) at t — and hence this equation give 
us correctly the time dependence of the electric field for all times. Finally 
having the electric field operator, we can obtain the magnetic field using 

ik x E_(k, t) = B_(k, t) and then the polarization densities employing the 
constitutive Eqs. fl33|) and f )38|) . 
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6 The magnetodielectric media with only one 
response equation 

In the previous sections we have assumed the magnetodielectric medium 
under consideration is one for which the distinction between polarization 
and magnetization physically is possible and the total current induced in 
medium can be separated into electric and magnetic parts. It is remarkable 
that the use of polarization and magnetization , separately, to describe the 
response of a medium to electromagnetic field is useful for media which are 
not spatially dispersive or which have a particular form of spatial dispersion 
[2H] . There are in general media for which it is not clear how the separation of 
the total current induced in medium is to be made into electric and magnetic 
parts. The disturbances, caused by the electromagnetic field, in such media 
are described completely in terms of only a response equation relating the 
total current induced in medium to the electric field [26J. The quantization 
method outlined in the previous sections can cover such cases applying the 
Lagrangian 



L(t) 



du l d r 



-X w ■ X w — -u 2 ^^ ■ X u 



+ 



+ 



d 3 r 



-e E 2 

2 2/i 



du I d A r I d 3 r'f l:j (u,r,r')E\f,t)Xi(r',t) 



(60) 



which is the same as the Lagrangian ([I])-© with the exceptional that, now 
the medium is modeled by a single set of three dimensional harmonic oscil- 
lators Xu- In this case the Euler- Lagrange equation (fT4j) takes the form 



SL 



5L 







i = 1,2,3 



8 [Ai(k,t) J J 5{A*(k,t)^ 
^oA(k, t) + [xqEq ik ip(k, t) — k x ( k x A(k, t)\ = fi P(k, t) 



where now 



P(k,t)= I du I d 6 pl(u,k,p) - X^(p,t) 



(61) 
(62) 
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can no longer be interpreted as the electric polarization, but is related to 

the total current induced in the medium by j(k,t) = P_(k,t). Following the 
canonical quantization similar to the previous sections we obtain the only 
response equation of the medium as 

P(k, t) = P N (k, t) + e jT" dt' j d 3 q x(k, q,\t\-t')- E(q, ±t') (63) 

instead of the two constitutive equations f )33|) and fj40|) . where P_ N (k,t) is 
given by ( 1371) and the response tensor 

±J o ™du^fd s pl(u,k,p)-f(u,q,p) t>0 
x (k,q,t)={ (64) 

t < 

is related to the conductivity tensor of the medium by |= = g_. This is one 
of the two alternative descriptions of the response of a magnetodielectric 
medium for which the distinction between polarization and magnetization is 
not possible [2E]. Finally following the same calculations as before sections 
one can obtain the space -time dependence of electromagnetic field operators 
in terms of the response tensor x an d the dynamical variables at t = 0. 
The calculations in this case are precisely the same as before sections except 
that the quantities Y^,]]^, M, x m ■> 9 are identically zero. 



7 Summary and conclusion 

By modeling an anisotropic polarizable and magnetizable medium with two 
continuous collections of three dimensional vector fields, we have presented a 
fully canonical quantization for electromagnetic field in the presence of such 
a medium. In this model the polarization fields of the medium did not enter 
directly in the Lagrangian of the total system as a part of the degrees of 
freedom of the medium, but the space-dependent oscillators modeling the 
medium solely was sufficient to describe both polarizability and the absorp- 
tion of the medium. The interaction of the medium with the electromagnetic 
field explicitly was given and both the Maxwell's laws and the constitutive 
equations of the medium were obtained as the consequences of the Heisen- 
berg equations of total system. The electric and magnetic polarization fields 
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of the medium could be deduced naturally in terms of the dynamical fields 
modeling the medium. Some space dependent real valued tensor coupling 
the medium with electromagnetic field were introduced. The coupling tensor 
had an important role in this quantization scheme, so that the electric and 
magnetic susceptibility tensors of the medium were determined in terms of 
the coupling tensor. Also the noise polarizations were expressed in terms 
of the coupling tensor and the dynamical variables describing the degrees 
of freedom of the medium at t — 0. In the free space and in the case of a 
non-absorbing medium the coupling tensor and the noise polarizations tend 
to zero and this quantization method is reduced to the usual quantization in 
these limiting cases. As a tool, we have used the reciprocal space to quan- 
tize the total system. It was shown the temporal backward and forward 
Laplace transforms of the electric field obey some algebraic equations with 
source terms that could be determined in terms of the dynamical variables 
of the system at t — 0. By introducing the Green function of these algebraic 
equations we were able to determine the Green function up to an arbitrary 
accuracy using a perturbation method which may be useful in some real cases 
with weak polarizability. The time-dependence of electric field was derived 
for both negative and positive times. Finally the model was modified to 
the case of a medium for which the distinction between the polarization and 
magnetization is not possible. 
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